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I. INTRODUCTION 



Much of our knowledge on the electronic properties of materials has been accumulated 
over the years by studying the spectra of electrons inelastically reflected from the surfaces 
of solid samples. The corresponding technique is usually referred to as (reflection) electron 
energy loss spectroscopy (EELS) [1]. By bombarding the solid sample by a monochromatic 
beam of electrons and measuring their energy loss and their deflection angle in the reflection 
mode, detailed information on the collective excitations of surfaces can be extracted, such 
as the excitation energies, the lifetimes, and dispersions. The energy loss of the impinging 
electron upon traversing the sample results in a variety of excitations including phonons and 
plasmon excitations, interband and intraband transitions, and inner-shell ionizations (the 
latter is particularly useful for detecting the elemental composition of a material). 

In the past two decades qualitative advances have been achieved in the capabilities of 
the electron-pair spectroscopy, also called the (e, 2e) spectroscopy, in which one studies the 
energy and the angular distributions of two electrons emitted simultaneously from a surface 
following the impact by one electron IJQ. In essence, the (e, 2e) spectroscopy investigates 
a specific process that, among others, also contributes to the EELS signal: the ejection of 
a secondary electron induced by the interaction of an impinging, primary electron with the 
surface. The basic information that one can obtain in this way is as follows: (i) the surface 

nn 

one-electron spin-resolved spectral function |J-|8[, (ii) the mechanisms of electron-electron 



collisions at surfaces 



14j |. and (iii) the surface dielectric function ISMlTj]. Concerning the 



spectroscopic data on the surface one-electron states, the (e, 2e) method is close to angle- 
resolved photoemission spectroscopy (ARPES) [isl, 19 1 because both measure the energies 
and the wave vectors of the surface electrons. An important difference between the two 
techniques is that they involve different surface one-electron transitions. This is due to the 
fact that in the (e, 2e) case the Coulombic force imposed by the projectile on the surface 



electrons is parallel to the momentum transfer [20|] while in the ARPES case the imposed 

electric force is, by definition, perpendicular to the momentum transfer, that is, to the 

photon momentum. Note in this context that, in contrast to ARPES, in the (e, 2e) method 

one can vary the momentum transfer in a wide range. Another marked feature is the high 

surface sensitivity of the (e, 2e) method, especially in the grazing-incidence mode, which 

I ' j I 

makes it very promising for exploring the Shockley [21|] and Tamm |22| electronic states 



found at the surfaces of various materials. 



The (e, 2e) measurements on LiF films deposited on a Si(OOl) surface [23|| gave the first 
evidence that the (e, 2e) spectroscopy is capable of providing insight into the secondary- 
electron ejection assisted by a collective excitation in the surface electronic system. A 
recent experimental study [15|] evidences a notable enhancement of the electron-induced 
secondary-electron emission from an Al(lOO) surface at the electron-energy- loss values that 
are equal to the bulk and the surface plasmon excitation energies. The above findings point 
out the potential for studying directly the plasmon-assisted electron-electron collisions and, 
in particular, the mechanisms of the plasmon decay at surfaces with the (e, 2e) method. 
Moreover, due to its remarkable surface sensitivity the (e, 2e) spectroscopy can be very 
successful in illuminating the properties both of conventional surface plasmons supported 
by metallic surfaces and of the low-energy acoustic two-dimensional plasmons that have 
been recently predicted 25, 26| and observed js^]. This calls for the relevant theoretical 
framework that consistently incorporates the effects of the surface dielectric response into 
the treatment of the (e, 2e) process, a task that hitherto remained outstanding and will be 
treated in this work. 

Two main mechanisms of the secondary-electron ejection from metallic surfaces following 



the electron impact are possible 28|: (i) due to the direct scattering between the incident 
(primary) and the valence-band (or the conduction-band) electrons and (ii) due to the decay 
of the bulk and the surface plasmons excited by the incident electrons. The bulk-plasmon 
decay into a single electron-hole pair is governed by the interband transitions, which in the 
case of long-wavelength plasmons are practically vertical in the the reduced- zone scheme. In 
a iellium model the interband transitions are absent and, for instance, within the random 
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phase approximation (RPA) [29| the decay of a plasmon into a single electron-hole pair 
can take place only when its momentum exceeds some critical value, namely, when the 
plasmon-dispersion curve merges with the electron-hole continuum. The simplest scenario 
for the direct electron-ejection mechanism assumes a single electron-electron interaction via 
the Coulomb potential screened by the surrounding medium. Dynamical screening effects 
can result in a resonant enhancement of the potential at the energy transfers corresponding 
to the excitation of collective modes, such as bulk and surface plasmons. This feature can 
markedly manifest itself as a large increase in the yield of the secondary electrons when the 
energy loss of a charged projectile is resonant with the plasmon energy (see, for instance. 



Ref. Since the contributions from both mechanisms exhibit a resonant behavior at 

the plasma energies, their separation in the (e, 2e) experiment is not straightforward. On 
the theoretical side, simple estimates show that the ratio of the (e, 2e) rate due to a direct 
scattering to that due to the plasmon decay behaves as oc F^^^, where Tpi is the plasmon line 
width. Therefore, the sharper the plasmon resonance is, the more dominant the role of the 
direct scattering over that of the plasmon decay is. 

In the present work, we consider theoretically the plasmon-assisted (e, 2e) collisions at the 
surfaces of the independent-electron metals and focus on the contribution from the direct- 
scattering mechanism. One of the key ingredients that determine the (e, 2e) rate in this case 
is the dielectric response of the metallic sample. The presence of the surface brings about 
a major complication as the response function undergoes a sudden change at the metal- 
vacuum interface. We address the problem of the surface dielectric response by considering 
two well-known approaches. The first one is the so-called specular- reflection model (SRM), 



which was first introduced in Ref. 



30j to study surface plasmons. The other is based on RPA 



with infinite surface barrier (RPA-IB) for electrons in metal 3l|. Below we incorporate both 
approaches in numerical calculations for (e, 2e) from Al and Be. The choice of Al and Be is 
motivated by the fact that among the free-electron metals they exhibit respectively a sharp 
and a wide plasmon resonance (in terms of the ratio between the plasmon line width and 
the plasma energy in these materials). It allows us to inspect and numerically illustrate the 
role of plasmon resonances on the dynamical screening of the electron-electron interaction. 

This paper is organized as follows. In Sec. [Ill we present basic formulas and approxi- 
mations for the considered process. In Sec. IIIIl models of surface dielectric response are 
discussed. Then, in Sec. IIV| numerical calculations for Al and Be are presented and ana- 
lyzed. Sec. |V] is devoted to a comparison of the present theoretical formulation with the 



recent experimental measurements on Al 15|, |l6|. Finally, conclusions are drawn in Sec. IVII 



Unless otherwise stated, atomic units (a.u., e = h = rrie = I) are used throughout. 



II. GENERAL FORMULATION AND BASIC APPROXIMATIONS 

We consider the process where, following the impact of a fast impinging electron with a 
momentum kg and energy Eq, two electrons are emitted from the surface of a semiinfinite 
solid with momenta k^, kg and energies Eg, E^. (see Fig. [1]). Hereafter the subscript s (e) 
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stands for the scattered (ejected) electron. All the energies are measured with respect to 
the vacuum level, so that in vacuum we have the dispersion 



kj = ^/2E, (j=0,s,e). 
The rate of the discussed reaction is determined by the so-called (spin-averaged) fully 



differential cross section (FDCS) [32 1 
da k^k, 



|^|(k,ke|T|ko«) + (kek,|r|koi)| 



dEsdVt.dEJVte {2TxYk 

+^|(k,k,|T|koz) - (k,k,|T|koz)p| 5{E, + E,-Eo- e,). (1) 

Here we specified the directions of the momenta of the emitted electrons by the solid angles 
Qs/e- The state vectors jk^ke) and |koi) describe, respectively, the two final-state electrons 
with asymptotic momenta k^, kg and the initial state consisting of the projectile state with 
momentum ko and the valence-band state The sum is taken over all occupied one- 
electron states of the surface with energy ei = Eg + E^. — Eq. The operator T is an effective 
transition operator that induces the (e, 2e) process and is assumed to be spin independent. 
In the frozen-core approximation it has the formal structure 

T = K + + (K + K + W)GUEtot)iVs + W), (2) 

where Vs,Ve, and W are effective (in general, optical) electron-solid and electron-electron po- 
tentials, respectively, and G^^{Etot) is the retarded two-electron propagator in the potential 
Vs + Ve + W at the total energy Etot = Es + Ee- The latter satisfies the Lippmann-Schwinger 
equation 

GUEtot) = GtiEtot) + G^{Etot){Vs + Ve + W)GUEtot), (3) 

with GQ{Etot) being the free two-electron propagator. 

In what follows, we treat Eq. ([2]) only to the first order in the electron-electron interaction 
W. Such a procedure is usually justified by the choice of the kinematics such that Eq ^ AE 
{AE = Eq — Eg) as well as by the screening of the electron-electron interaction due to 
the surrounding medium. This then amounts to the distorted wave Born approximation 
(DWBA) 

T = [1 + (K + V,)GtMot)] W[l + Gt,AEtotm (4) 



where Gf^{Etot) is the two-electron propagator in the potential Vs + Ve and is given by the 
solution of the following Lippmann-Schwinger equation: 

Gt,{Etot) = G+{Etot) + G+{Etot){Vs + Ve)Gt,{Etot). (5) 



Taking into account Eq. (jlj), we can present FDCS ([I]) as 



da k^kf. 



where 



degdVLgdeedVLe (27r)^A;o 1 4 

x5{Es + E,- Eo-ei), (6) 

xLtV) = ^'''^ + / dr'e'^°^'v{v')g+{v', r; Eq), (7) 

xL;^(r) = e^"^^'^ + j dr'e'^^^'v{r')g-iT', r; i?,), (8) 

with j = s,e and V^- = v, and (7^ ((7^) being the retarded (advanced) one-electron Green's 
function in the potential v. In the case of solids with a translational symmetry parallel 
to the surface, the states Y y^\ can be calculated within the dynamical low-energy electron 



diffraction (LEED) theory |33l|34|. 



Employing the surface dielectric function e for the description of the effect of the screening 
on the bare electron-electron interaction Vee, we get 

W{r, r'; u) = j dr"e-\r, r"; uM", r'), v^r", r') = ^1-^, (9) 

where is the inverse dielectric function, and the energy argument CO (cij — Eq — Eg or 
CO = Eq — E^, depending on the final state of the incident electron) accounts for the dynamical 
screening effects. 



III. SURFACE DIELECTRIC RESPONSE 

For our purposes we need the inverse dielectric function that is derived from the 
dielectric function e according to 

J dr"e{r, r"; io)e-\r", r'; ^) = J dr"e-\r, r"; uj)e{r", r'; uj) = S{r - r'). (10) 



Let us assume that the z axis is perpendicular to the surface and that the sohd sample fills 
the 2; < region (see Fig. [1]). If the sample is crystalline, then 

e(r,r';w) = e(r + A,r' + A;w), ^-^(r, r'; w) = ^-^(r + A, r' + A; w), (11) 

where A is the lattice vector parallel to the surface. Thus, e and can be presented as 

.(r,r»=5^e^(«^-«'^') / ^ .o,g'(Q, ^, ^»e^^(^-^'\ (12) 

J (2vr)^ 

l^tBZ 

.-(r,r';c.) = e^^^^'^'^'^ J ^ .-,,(Q, z»e^Q(^-^'), (13) 



G,G -^,tg2 



where r = (R, -z) and r' = (R',2;'), G and G' are the surface reciprocal lattice vectors, and 
the Q integration is carried over the first Brillouin zone (1*^* BZ) of the surface reciprocal 
lattice. Using Eq. f|T3|) . we get for the screened Coulomb potential Qj 

W{r, ro; c) = e^(«^-«'^o) J ^ J dz'e-^]M ^ -)eHQ+^'ll^'--l. 

(14) 

If we neglect the crystalline effects on e and parallel to the surface, then 

£g!g'(Q' ^) = '^G,G'£^~^(Q + G, z'; w), (15) 

where it is assumed that Q 6 BZ, and hence 

y dz'e-\Q,z,z'-u;)e-'^\^'-^^\. (16) 

Representations (IT^ and ffT6|) can be particularly useful, when employing the following 
expansion of the electron states in Eq. ([6]): 

x'^\r) = J2ci^\G,z)e^^^^^^^, (17) 

G 

where k = (K, k^). 

If the sample is a free-electron metal (for instance, Al or Be), then the model of a 
degenerate electron gas where electrons move in a positive ionic background and are bounded 
by a surface potential barrier is commonly applicable to mimic its dielectric properties. 
Below, we briefiy sketch two possible approaches for describing the dielectric response of 
such jellium-like systems. 



A. Specular-reflection model 



The problems involving the interaction of charged particles with plane-bounded solids 
are often treated using SRM. This model assumes the surface barrier to be impermeable for 
the electrons in the solid, so that they are specularly reflected at the surface. The reflection 
process is described in a classical fashion, in particular the interference between the incoming 
and outgoing components is neglected. Within SRM the potential created by the external 
charge distribution p(r, t) in the vicinity of a surface is given by (see, for instance, Ref. jssj 
and references therein) 

V{r,t) = V4r,t)e{z) + V4r,m-z), T4(r,t) = J J ^V^{Q, z;u;)e^^^''--'\ 

(18) 

where 

V±{Q, z- u) = 47r[[/±(Q, z; u) T p(Q, a;)i^±(Q, z; u)], (19) 
with [using the notation q = (Q, g^)] 

[/.(Q,.;^) = /^^l^i^le-, £/_(Q.^;-) = / pl^^e'^.^, (20) 
J 211 J 27r g^£:6(q,w) 



p±(q,aj) = J dr J dt p{Il,±\z\,t)e~ 



■i{qr— ojt) 
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MQ,z;u)= [ ^^e'^- = e-«l^l, U-{Q,z-u;)= [ ^ , f , e^^--, (21) 
J 211 J 2% g^£:fe(q,a;) 

i + EsKSl-,'^) 

In Eqs. ( I20l) and ( 12T1) . ^^(q, w) is the bulk dielectric function, i.e., that of an infinite 3D 
system.^ The quantity SgiQl^uj) occurring in Eq. (122|) is the so-called surface dielectric func- 
tion [31i] . 

To utilize the SRM approach in the present quantum- mechanical treatment, the ex- 
ternal charge density p(r, t) must be replaced by the corresponding operator p(r, t) = 
e~^^°^p{r)e^^°^, where Hq is the Hamiltonian associated with the incident electron. Clearly, in 
this way the frequency u equals the energy transfer AE = Eq — Es^e)- Since p(r) = 5(r — fq), 
where rg is the position of the incident electron, it is straightforward to show that the 
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resulting expression for the inverse dielectric function in Eq. ( !T6|) is 



e-\Q, z, z'- u) = e{z) <! 6{z - \z'\)e{zo) - , M ^'Mzo) - «:(Q, z'; a;)0(-^o)] 



+ei-z) I [fi:(Q, z + z'; cu) + k{Q, z - z'; uj)]e{-Zo) + j 



2z/_(Q, z; u) 



x[S{z')e{zo) - K{Q,z';io)e{-zo)] }, (23) 



where 



k{Q,z;u)= I (24) 



Note that the inverse dielectric function ( l23l) depends on whether the incoming electron is 
inside {zq < 0) or outside (2:0 > 0) the solid. 



B. Random phase approximation 

RPA constitutes a reasonable framework for describing the dielectric response of a degen- 
erate electron gas. The bulk dielectric function was first derived within this metho d by 
Lindhard 36| . For the case of a semi- infinite geometry a very useful study is from Newns 3l| , 
who considered the dielectric response of a semi-infinite ideal metal within RPA assuming 
an infinite surface barrier (RPA-IB), which is in the spirit of SRM. Using the approximation 
of a specular electron reflection, Bechstedt et al. js?}] derived expressions for thj^^screened 
Coulomb potential and the inverse dielectric function. However, Horing et al. [38] noticed 



that the result of Ref. [37| is incapable of describing correctly the image field as a part of 
the dynamically screened Coulomb potential. In their calculation, Horing et al. utilized the 
potential solutions obtained by Newns |31| within the RPA-IB model. Employing the math- 
ematical model delta-function potential unscreened interaction and neglecting the 
nondiagonal elements in the density-response matrix, they calculated the inverse dielectric 



function as 



e-\Q, z, z'- 00) = e{z) !^6{z - z') - ^-^^—^[6{z') - 2k{Q, z'; oo)e{-z')]j 

+Q{-z) I [Kid, z + z'- u) + «:(Q, z - uj)]Q{-z') + '^-}^]'''^\ 
x[6{z') -2k{Q, z';uj)Q{-z')]^, (25) 

where the functions k(Q, z; u), z^_(Q, z; u), and ^^(Q, oj) are the same as in Sec. IIII A[ 

Expressions fl23p and fl25p have a similar structure. Moreover, one can formally de- 
duce Eq. (!25l) using the SRM approach in the case of the model delta-function potential 
Vee{r,ro) = 6{r — tq). However, the RPA-IB result fl25l) is supposed to be applicable in the 
case of Coulomb-like potentials as well (see Ref. j38|] for detail). The latter feature makes 
the two approaches nonequivalent. 



IV. RESULTS AND DISCUSSION 

In this section we present and analyze the numerical results for the correlated electron-pair 
emission from Al and Be surfaces. To be specific, we consider the kinematics of the recent 

energy loss was measured at 

the incident energy Eq = 100 eV and the angle 6q = 30° in the specular reflection mode, that 
is, 9s = 30°, in coincidence with the secondary electron ejected at the angle 9e = 60° (see 
Fig. [1]). We neglect the crystalline effects in our calculations and construct the one-electron 
states in Eq. (E]) in the context of the jellium model. The details of this procedure are given 
in the Appendix. The surface dielectric response is taken into account within the SRM and 
RPA-IB approaches described in the previous section. Each of them depends on the model 
of the bulk dielectric response. Below we outline two approximations for the bulk dielectric 
function of a free-electron metal that were employed in the present calculations: (i) the 



Thomas-Fermi (TF) |39| and (ii) the hydrodynamic (HA) [40|| approximations. In spite of 
their relative simplicity, they efficiently mimic the basic features pertinent to the static and 
dynamical bulk screening effects in metals. In addition, their use in the calculations makes 
the numerical implementation more transparent and controllable. 
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1. Thomas-Fermi approximation 



This well-known model neglects the dynamical screening effects and assumes the dielectric 
function to be of the form 

= 1 + — , (26) 
where A is the screening constant. For a free-electron isotropic gas one has 



vr 

where ep and kp are the Fermi energy and momentum, respectively. Using the Wigner-Seitz 



radius r^, we have 



(r, = 2.07 for Al and r, = 1.87 for Be). 

Within the TF model the functions in Eqs. (1231) and (1251) are given by 



«:(Q, ^; uj) = 6{z) - ^ e'^^^ u.{Q, z;u;) = ^ e'^l^ ^.(Q, oo) = |, (27) 



where A = ^/Q^ + X^. 



2. Hydrodynamic approximation 



Since its introduction almost 80 years ago 40[, the hydrodynamic approximation has 
proved to be very useful in describing the electrical transport and the optical properties of 
conductors. The main advantage of this model lies in the simplicity of accounting for the 
spatial dispersion, 

Sbiq, uj) = l+ ^ (28) 
where the frequency of the bulk plasmon mode Ub is given by 



Ub = yi-nn = w — , 
y 

with n being the electron density. The parameter (3 depends on cu, such that 

T 



-Vf, u <^u 
/3=<! ^1 (29) 
-Vf, w > z/. 
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The above low- and high-frequency hmits can be reproduced with the "interpolation for- 



mula" 



4l| 



^„ fo; -I- ^iu r, , , 

p2 ^ 5 3_^2 ^ ^3Q) 

CO + iu 

The collision frequency u can be estimated as ~ Tpi (Tpi = 0.53 eV for Al and Tpi = 4.7eV 
for Be). Note that the TF result (126|) derives from Eq. fl28l) in the limit — )■ 0. 
Within the HA model the functions in Eqs. ( !23l) and ( |25|) are given by 



k(Q, z; cu) = 6{z) — 



ujI e-^l^l 



i/_(Q, 2;; w) 



2/32 A ' 



+ iu) — ojI a ci;(ci; + iu) — ' 



where A = —i/S ^ ^yuJ{uJ + iu) — ujI — /S'^Q'^ (here the first branch of the square root function 
is assumed, which yields Re(A) > 0). 



A. Aluminum 



Fig. [2] shows the correlated electron energy distribution of the scattered and the ejected 
electrons calculated using the SRM surface dielectric function (1231) . Marked differences 
between the TF and HA models can be seen. The maximum of the intensity in the TF case 
is located at small values of -Eg ~ 1 eV when the energy loss AE = Eq — Es exceeds the work 
function $ = 4.3 eV by approximately the same amount, which indicates that the ejected 
electron originates from the initial state close to the Fermi level ep- With increasing AE 
the intensity decreases and a tendency can be observed: at a given energy-loss value the 
electrons are preferably ejected with energies close to the threshold E^ = AE — $, which 
corresponds to ei = ep. The latter is readily explained by the fact that the density of states 
at the Fermi level is maximal. In contrast to the TF model, in the HA case the intensity 
is strongly peaked around AE ^ lleV which is slightly above the surface plasmon energy 
Us = cob/V2 = 10.5 eV. The reason for such a resonant behavior is due to the poles of the 
function (1251) when 1 + Es = 0. If neglecting the plasmon dispersion and the damping effects 
in the HA model fl2S]) . that is, /3 = and 1^ = 0, one finds that, according to Eq. ( 1HT]) . the 
pole is located exactly at AE = Us- The plasmon dispersion and the damping effects are 
responsible for the shift and for the finite width of the observed resonant peak. 
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The numerical results for the correlated electron energy distribution using the RPA- 
IB surface dielectric function f l25|) are presented in Fig. |3l The TF distribution is very 
close to that when using SRM. However, in the HA case the distribution differs from the 
analogous one in Fig. |5J Namely, in addition to the peak associated with the surface plasmon 
there appears another, more pronounced feature at approximately AE = 15 eV. Taking into 
account that Ub = 14.9 eV, it is clear that the feature is due to the dynamical screening 
effects related to the bulk-plasmon mode. The absence of the bulk-plasmon peak in the 
SRM case and its appearance in the RPA-IB case follows from the comparison of Eqs. ( 123|) 
and fl2Sl) . The bulk screening effects in these models are associated with function (^^. 
Thus, within the SRM (1231) it comes into play when the incident electron penetrates inside 
the metal {zq < 0), while within the RPA-IB f l2^ it becomes already relevant when the 
incident electron is still moving in the vacuum. A short inelastic mean free path for the 
incoming electron results in small penetration lengths, thus strongly restricting the bulk 
contribution in the SRM case. 

Fig. m compares the so-called five-fold differential cross sections (5DCS), which derive 
from FDCS ([6]) upon integrating over the ejected electron energy Eg, using the SRM and 
RPA-IB approaches. Clearly, 5DCS characterizes the dependence of the ejected-electron 
yield on the energy loss AE in the considered collision geometry. In accordance with the TF 
results in Figs. |2]and[3l the TF ejected-electron yields in Fig. HJare very close to each other, 
starting to grow at the threshold value AE = $, exhibiting maximum at AE = 7 — 8eV, 
and then smoothly decreasing down to at AE ^ 16 eV, which is due to the restriction 
|Ks + Ke — KqI < kp imposed by the conservation of the electron-pair momentum parallel to 
the surface. The HA results are orders of magnitude larger than the TF ones, indicating the 
strength of the dynamical screening effects. In accordance with Figs. |2] and |3l the HA results 
in Fig. m using SRM exhibit only one peak associated with the surface-plasmon mode while 
those using RPA-IB two peaks which are due to both the surface- and the bulk-plasmon 
modes. The intensity of the surface-plasmon peak in the SRM case exceeds that in the 
RPA-IB case by almost two orders of magnitude. The position and the intensity of the 
bulk-plasmon peak in the RPA-IB case is mainly determined by the interplay between the 
plasmon pole in the HA bulk dielectric function and the kinematical effects related to the 
conservation of the electron-pair energy and the surface-parallel momentum. 
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B. Beryllium 



Figs. |5l O and [7] present numerical calculations for the same situations as in Figs. |2l |3l 
andlU respectively, but for Be. In the HA case, the differences compared to the results for Al 
can be attributed to a much wider plasmon resonance in Be. In particular, we find that the 
footprints of the plasmon modes in the HA results are much weaker and broader. The role of 
the plasmon linewidth Tpi is clearly seen when comparing the HA results for Al and Be using 
SRM. The FWHM of the peaks observed in the HA results within SRM in Figs. H and [7] 
are given by approximately Tpi{Al) = 0.53 eV and rp;(Be) = 4.7eV, respectively. Moreover, 
when using RPA-IB, the plasmon features can hardly be identified in FDCS (Fig. [H]). And 
they manifest themselves for some cases only in 5DCS (Fig. [7]). In contrast, the TF results 
for Be are rather close to those for Al, particularly in magnitude. Though the differences 
between the TF and HA results in Figs. [5] and [6] are not as large as in the case of Al, the 
effects of dynamical screening are still quite strong. This conclusion follows from the fact 
that the HA results for 5DCS are about four orders of magnitude lager than the TF ones 
(see Fig. ED. 



V. THEORY AND EXPERIMENT 

Here we discuss how the predictions o, 
results of the recent experiments on Al 



le present theoretical approach compare to the 
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16| . The kinematica! 



regime of Ref. 



15| was 



specified in the previous section. The measurements in Ref. [l6| were carried out for a 
normal incidence {6q = 0°) of the projectile electron with an impact energy of 500 eV. 
The scattered and the ejected electrons were detected at the polar angles Og = 6e = 60° 
(see Fig. [1]). In both studies, the secondary-electron spectra were recorded in coincidence 
with the primary electron having undergone an energy loss that is equal to the bulk- or the 
surface-plasmon frequencies, i.e., when AE = Ub and AE = uog- In coincident measurements 
of the secondary-electron spectra in Ref. [3], further data were reported for the energy-loss 
values of 25, 30, 40, 45, and 150 eV. 

For a detailed quantitative comparison of our theoretical results with the discussed mea- 
surements one should consider the following aspects. First, the experimental data are re- 
ported on an arbitrary intensity scale, which means that the absolute values of the FDCS are 
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not determined. Second, one must bear in mind that the experimental data are broadened 

n 

by a finite energy and angular resolution whose effective values in Ref. [15| were given to be 
1 and 1.2 eV for the primary (scattered) and the secondary (ejected) electrons, respectively, 
while in Ref. jl6| the overall energy resolution was reported to be about 5eV. Below we take 
into account the energy-broadening effect by convoluting our theoretical calculations with a 
Gaussian energy distribution. 

Figs . [8] and [9] show the results of the numerical calculations for the FDCS in the kinematics 
of Ref. [l5| , which were convoluted with the 2D Gaussian function 



P(S:, K) = exp 



2(7? 2ai 



where FWHM<, = leV and FWHMe = 1.2 eV. As remarked above, the absolute values 
were not measured in the experiment, but the scale for the two coincidence spectra (at 



AE = uJb and AE = Ug) in Ref. [15| is the same. In order to place the experimental results 
on a common intensity scale with the theory in each panel of Figs. |8] and [9l they are nor- 
malized in such a way that the maximal experimental and theoretical values of the FDCS 
corresponding to the AE = Us case are the same. In terms of the positions of the peaks 
in the secondary-electron spectra, all calculations agree reasonably well with experiments. 
However, in contrast to the experimental results, the theoretical spectra exhibit apprecia- 
ble intensities in the region E^. < 4eV. With regard to a comparison of the intensities in 
the AE = Ub and AE = Ug cases, only the RPA-IB model using the HA bulk dielectric 
function provides a reasonable agreement with experiment. Indeed, the TF approximation, 
both within RPA-IB and within SRM, predicts the AE = Us results that are by an order of 
magnitude larger than the AE = Ub ones, while the SRM model using HA yields an even 
much larger discrepancy (almost three orders of magnitude). Thus, it can be concluded that 
within the present theoretical approach the best overall agreement with the experimental 



data of Ref. |15| is found when using the RPA-IB model that involves the HA bulk dielec- 
tric function. This comparison hints at the suitability of the measurements to asses the 
reliability of the employed dielectric response. We note, however, that further comparisons 
with measurements at different geometries, as well as on different samples, are necessary for 
conclusive statements on the quality of the discussed approximations. 



The experimental conditions chosen in Ref. 16| are beyond the validity range and the 



initial scope of our theoretical treatment. In fact, in this case the theoretical FDCS vanishes 
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for all AE that were chosen in the experiment. The reason for this is that in the theory we do 
not account for coherent or incoherent multi-scattering events that randomize the momenta. 
As a result we are bound to the kinematical limitation |Ks + Kg — Ko| < kp that stems 
from the assumption of an electron ejection upon a direct electron-electron scattering. In 
principle, one may lift this limitation within the present theoretical approach by taking into 
account the effects of electron-pair diffraction 
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23( 1 and surface roughness. In our opinion 



this would be, however, insufficient to explain the dominant, broad peak- like structures 
with a falling edge at about Ee = Ub — ^ that are observed in the measured coincidence 
spectra when AE exceeds largely the Uh value. Clearly, due to the energy balance, these 
features can not be accounted for within the picture of only one inelastic electron-electron 
collision involving either plasmon decay or dynamical screening effects. The experimental 
findings [16] thus hint at the existence of multiple inelastic scattering in the kinematics under 
study. In this context we remark that the influence of these multiple-scattering processes 
on the pair correlation functions was discussed and experimentally verified for a LiF sample 
in Ref. In Ref. [isl it was argued that multiple-scattering processes can be viewed as 
a Markov chain, with independent successive inelastic collisions leading to the excitation of 
the surface and the bulk plasmons followed by their decay into single electron-hole pairs (see 
Ref. [16] for detail). We add that one cannot rule out the scenario in which the individual 
inelastic events in the Markov chain take place due to direct electron-electron scattering 
resonantly enhanced at the energy transfers close to the surface- and the bulk-plasmon 
frequencies. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we considered theoretically the electron energy-loss process accompanied by 
the ejection of the secondary electron from a clean metallic surface. Our study was focused 
on the effects of the surface dielectric properties in the discussed process. Restricting the 
analysis by the mechanism of a direct electron-electron scattering, we took into account the 
modification of the bare Coulomb potential due to the surrounding medium by means of the 
inverse dielectric function. Using two models of the inverse dielectric function of the surface, 
SRM and RPA-IB, we performed numerical calculations for Al and Be with and without 
accounting for the dynamical screening. We found that the results exhibit clear plasmon 
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features in the case of Al, while in the case of Be these features are much less pronounced, 
which can be attributed to a wider plasmon resonance in Be. For both metals, very strong 
dynamical-screening effects were found. 

In the present theoretical analysis we did not address the issue of the plasmon-decay 
mechanism which also contributes to the considered process. Studying this specific mecha- 
nism with the (e, 2e) method can offer an opportunity for the detailed investigation of the 
plasmon-decay channel into a single electron-hole pair. For instance, the recent experimen- 
tal results from Al [15 showed clear surface- and bulk-plasmon features in FDCS that were 
interpreted in Ref. [l5j within the plasmon-decay model. However, as shown by our RPA- 
IB HA calculations in Fig. |8l these features can be interpreted within the direct-scattering 
model as well. As stated in the Introduction the role of the plasmon-decay mechanism is 
expected to be more significant in free-electron-type metals exhibiting a broad plasmon res- 
onance, such as Be. However, to our knowledge, no coincidence (e, 2e) measurements on 
such materials have been published so far. 

Thus, further experimental studies on the plasmon-assisted (e, 2e) collisions at metallic 
surfaces are desirable in order to shed more light on the roles of the discussed mechanisms. 
On the other hand, to account for experimental observation theory must go beyond the 
framework of a single inelastic collision. 



ACKNOWLEDGMENTS 

We are grateful to Jiirgen Kirschner, Arthur Ernst and Sergey Samarin for useful dis- 
cussions and to Giovanni Stefani and Alessandro Ruocco for valuable comments on their 
experimental studies. This work was supported by DAAD and SFB762. K.A.K. also ac- 
knowledges support from RFBR (Grant No. 11-01-00523-a). 



Appendix A: One-electron states 

Evaluation of the FDCS (jS]) requires the knowledge of the one-electron states that cor- 
respond to the incoming (incident), bound, and outgoing (scattered and ejected) electrons. 
Let us consider a semiinfinite metallic sample filling the space in the negative z direction. 

17 



Within the jelhum model, the effective one-electron potential V is a steplike one, i.e.. 



V = -Voei-z). (Al) 

For a clean metallic material one has 

Vo = eF + $, (A2) 

where $ is the work function. The wave function and the energy of the electron bound in 
the metal are thus given by 

Xk(r) = e^^^ {B{h)e-^'e{z) + [e^^^^ + A{K)e-'''^^]e{-z)} , = ^(K^ + k^) - Vo, 

(A3) 

where 

^(''^-"Brr «(*')-^- ^'-V2v^. (A4) 



and y/2Vo >h>0. 

Within the SRM and RPA-IB approaches the surface barrier is impermeable for electrons 
in the metal. This feature can be taken into account by calculating the wave function in 
Eq. (1A3P under the assumption Vq = oo, which yields 

Xk(r) = 2e^^^sin(M0(--2)- (A5) 

The incoming and the outgoing electron states with momentum k = (K, k^) and the 
energy = (K^ + k'^)/2 in the potential (lAip are given by 

xi;\r) = e'^^ {[e'"^^ + R{k,)e-''^']Q{z) + D(fc,)e''=^^e"^e(-z)} {k, < 0), (A6) 
xL^V) = e""^^ {e''^'Q{z) + [Ai(fc,)e*'=^^ + A2(fc,)e-*'=^^e"^e(-z)} {k, > 0), (A7) 



where 



kz + k'^ — ia kz + k'^ — ia' 

, _ k'^ + kz + ta A (u \ _ K- kz- 
~ 2k' ' ^2lfc^j - : 



K\ 



(A8) 
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with Voi being introduced as an imaginary component of the potential, 

V = {-Vo + tVo,)&{~z), 

to mimic the damping of the electron waves inside the metal. Using the concept of the 
inelastic mean free path A, it can be estimated as 

k' 



One can calculate A, for instance, from the parametrization formula by Seah and Dench 42| . 
which in the case of the elemental materials reads 

A= (^^^^^, +0-4Va3(E, + $), (A9) 

where the electron energy E'k + $ is measured in eV. The average thickness of the monolayer 
a measured in nanometers is given by 

= AlO 

pNa 

where M is the atomic weight, p is the bulk density (in kg/m^), and is Avogadro's 
number. 
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FIG. 1. (Color online) Schematic drawing of the electron-induced electron-pair emission from 
surfaces. 
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FIG. 2. (Color online) Correlated electron energy distribution from Al within SRM. The top panel 
corresponds to the TF bulk dielectric function, while the bottom corresponds to HA. The white 
solid line marks the energy threshold for the electron-pair emission, Eq — Eg — > where $ is 
the work function. 
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FIG. 4. (Color online) Dependence of the ejected-electron yields from Al on the energy loss within 
SRM and RPA-IB. The top (bottom) panel presents the results using the TF (HA) bulk dielectric 
function. 
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FIG. 7. (Color online) The same as in Fig. HI but in the case of Be. 
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FIG. 8. (Color online) Spectra of secondary electrons at energy losses equal to the surface- and 
bulk-plasmon frequencies. The calculations were carried out within the RPA-IB model using the 
HA (top panel) and TF (bottom panel) bulk dielectric functions. The experimental data points 



are borrowed from Ref. 
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FIG. 9. (Color online) The same as in Fig. [8l but within the SRM model. 
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